Renormalized waves and discrete breathers in /3-FPU chains 
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We demonstrate via numerical simulation that in the strongly nonlinear limit, the /3-FPU system 
in thermal equilibrium behaves surprisingly like weakly nonlinear waves in properly renormalized 
normal variables. This arises because the collective effect of strongly nonlinear interactions ef- 
fectively renormalizes linear dispersion frequency and leads to effectively weak interaction among 
these renormalized waves. Furthermore, we show that the dynamical scenario for thermalized f3- 
FPU chains is spatially highly localized discrete breathers riding chaotically on spatially extended, 
renormalized waves. 



The Fermi-Pasta-Ulam (FPU) lattice was introduced 
in their classical work [1] to address fundamental issues of 
statistical physics such as equipartition of energy, ergod- 
icity. The attempt to resolve the mystery of the FPU re- 
currence (the system did not thermalize as was expected 
but rather kept returning to the initial state [1]) has 
spurred many great mathematical and physical discov- 
eries, such as the celebrated KAM theorem and soliton 
physics [2]. Despite this remarkable progress, there are 
still fundamental open questions that are under vigorous 
debate [3], such as what is the route to thermalization 
and how to fully characterize thermalized /3-FPU sys- 
tem. Furthermore, in the last decade discrete breathers 
(DB) as spatially localized, time periodic lattice excita- 
tions were discovered [4]. Arising from energy localiza- 
tion in nonlinear lattices, they play important roles in 
many dynamics in fiber optics, condensed matter physics 
and molecular biology [5]. The existence of DBs has 
been addressed rigorously [6]. Important conceptual is- 
sues naturally arise, such as what is the role of DBs on 
the route to equilibrium [7] and how do they manifest in 
thermalization of the FPU system? Resolution of these 
issues will certainly provide deep insight into the funda- 
mental understanding of route to thermalization for gen- 
eral nonlinear physical systems. Most of the results re- 
garding DBs in /3-FPU chains have so far only addressed 
their behavior in the transient state of weakly nonlinear 
regimes before thermalization occurs [8,9]. 

In this Letter, we investigate the FPU dynamics in 
the strongly nonlinear limit. We demonstrate that, quite 
surprisingly, even for strong nonlinearity, the /3-FPU sys- 
tem in thermal equilibrium behaves like weakly nonlin- 
ear waves in properly chosen variables. Such behavior 
results from the collective effect of strongly nonlinear in- 
teractions effectively renormalizing linear dispersion rela- 
tion. This observation enables us to use a well-developed 
weak turbulence (WT) formalism [10] for the descrip- 
tion of the /3-FPU chains even in a strongly nonlinear 
regime. Furthermore, in addition to the nonlinear waves, 
we observe the DB excitations in the thermalized state 
of /3-FPU chains. Previously such DBs were observed 
only during transient stages towards thermalization [8]. 



Here we show via numerical simulation that DBs actually 
persist and coexist with renormalized waves in the ther- 
malized state. Thus, in the thermalized /3-FPU, there 
are two kinds of quasi-particle excitations, one localized 
in fc-space as renormalized nonlinear waves/phonons, and 
the other, localized in x-space as DBs. 

The /3-FPU chain is described by the Hamiltonian, 

H = H 2 +H 4 , (1) 
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H2 = ^ ^Pi + - *+i) 2 ; Hi = - ^(qi - qi+i) 4 , 

i=l i=l 

with periodic boundary conditions, where (3 is a nonlin- 
ear parameter, pi and qi , are the i-th particle momentum 
and displacement from the equilibrium position, respec- 
tively. In terms of Pk and Qk, the Fourier transforms of 
Pi and g;, the Hamiltonian becomes 

1 N ~ 1 

H = 2E W 2 +<&M] + V (Q)> ( 2 ) 
fe=i 

where ujk — 2 sin 2j£ is the linear dispersion frequency 
and V(Q) is the Fourier transform of H^. V(Q) is a lin- 
ear combination of various quartic products of Qfc's and 
Q£'s. We numerically integrate the canonical equations 
of motion for Hamiltonian (1) to study possible dynami- 
cal scenario of the FPU [11]. Since the thermalized state 
is our primary interest we chose random initial data [11]. 
Note that the behavior of the /3-FPU in equilibrium is 
fully characterized by only one parameter (3H [12] where 
H is the total energy of the system. We fixed H = 200 
and N = 128 (except for Fig. 6, where N=1024) and 
varied /3. 

First, for the FPU system in thermal equilibrium [13], 
with moderate and strong nonlinearities /3 = 1,8, and 32, 
we measured the power spectrum ( | cifc | 2 ) as a function of 
u>k, where aj~ = (Pk — iuJkQk) I sf^k an d (■) denotes the 
time averaging. Although it was expected that, for weak 
nonlinearity, (|a/c| 2 ) = T/ujk, where T is an effective tem- 
perature (Rayleigh- Jeans distribution for waves [10]), it 
is surprising to find that the same scaling holds even for 
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strong nonlinearities as shown in Fig. 1. To understand 
why (|afc| 2 ) scales as uj^ 1 even in the strongly nonlin- 
ear limit, we computed the u-k spectrum of dfc(i) as 
shown in Fig. 2(a). In the weakly nonlinear limit the 
spectrum would have resonant peaks along a curve given 
by the linear dispersion relation u) k = 2sm(kir/N). We 
observe that, when the nonlinearity is no longer small, 
the resonances move to higher frequencies, as the dis- 
persion relation is renormalizcd by the nonlinear part 
V(Q). The renormalizcd dispersion relations is indi- 
cated by the sine-like structure of the resonance peaks 
(the solid line in Fig. 2(a)). We verified numerically that 
the main contribution of the nonlinear potential energy 
( ~ 80%) comes from the terms QkQiQmQt constrained 
on k + I = to + s. We note that when k = to and I = s 
or k = s and I — to, these terms can be combined with 
the quadratic part of (2) to effectively renormalize the 
frequency u) k - More specifically, Hamiltonian (2) can be 
rewritten as the sum of a new renormalizcd quadratic 
part and a remaining nonlinear part: H = + H 4 
where H 2 = | E£L~i (\ p k | 2 + ^l\Qk\ 2 ) with the renor- 
malizcd dispersion relation, 
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Note that the frequency renormalization factor r\ does 
not depend on k. Fig. 2(b) shows how the renormalizcd 
frequency oj k depends on [3. The upper curve was pro- 
duced from the numerical spectrum the abscissas of 
the peaks of the sine-like curves were measured (e.g., from 



Fig. 2(a) for (3=l,u k ._ 



N/2 : 



=3.3 for N = 128). The curve 



with pluses was obtained using Eq. (3). Both seem to 
possess a scaling ~ /3 2 dependence (see Fig. 2(b)). The 
discrepancy between the analytical and numerical curves 
can be attributed to the fact that in the derivation of 
formula (3) we took into account only 4-wave processes. 
In principle, we can take into account higher order pro- 
cesses to obtain a more accurate analytical estimate of 
the renormalizing factor. 

To further study the renormalization of interactions 
we measured the ratio of the quartic to the quadratic 
parts of the energy before and after renormalization 
procedure (i.e., H4/H2 vs H4/H2) for different values 
of (3 with energy fixed (Fig. 2(c)). This figure also 
shows that the effective renormalized linear part becomes 
more dominant. Therefore even for strongly nonlinear 
regimes, with (3 as large as 128, the renormalized waves 
(or purely nonlinear phonons) have weakened interac- 
tions. Note that the resonance of a k (co) has a finite 
width (shown with the dashed lines in Fig. 2(a) which 
are the level of / \ak\ 2 (w)dw/ max TO \a k \ 2 for each k). We 
note that for k = N/2 (the highest mode in the sys- 
tem) the resonances become the broadest. According 
to WT [10], the energy exchange among waves occurs 
on the resonance manifold given by k\ + k 2 = fc 3 + k 4 



and Wfej + LUk 2 — w fe 3 + w fc 4 - Although one can show 
that there are no exact resonances in /3-FPU chains on 
the discrete lattice, the nonlinearity induced near res- 
onance interactions (i.e. \ujk 1 + uJk 2 ~ — w fc 4 | < $ 
where S is a resonance width) can occur. This allows us 
to use the WT theory if the interaction is considered to 
be weak. In order to characterize the system as weakly 
nonlinear waves, the renormalized waves are described 
by the new normal variables ak = (Pk — ioJkQk)/V^k- 
The relationship between bare ak and renormalized a k 
is a k = [{y/fj + l/y/fj)a k + {VV- Vv^-k] A Under 
the random phases approximation [14] for a k we have 
(|afc| 2 ) = 277/(1 + r/ 2 )(|afc| 2 )- Therefore the power spec- 
trum has scaling ui^ 1 for both ak and Sfe. (Note that 
(|afc| 2 ) is shown in Fig. 1). As a consequence, the effective 
temperature, T, for the renormalized waves, is related to 
the bare temperature T by T = [2ry 2 /(l + ?7 2 )]T. 

The renormalization picture is further corroborated by 
the time evolution of each individual wave ak- For the 
bare waves ak, there are large temporal modulations in 
both modulus and phase as a result of strong nonlinear 
interaction among these modes on the linear dispersion 
time scale (as shown in Fig. 3(a) and (b)). In weakly 
nonlinear systems, wave amplitudes and corresponding 
phases are expected to evolve slowly on the linear disper- 
sion time scale. For renormalized waves (Fig. 3(c) and 
(d))these modes indeed have characteristics of weakly in- 
teracting waves, with a small modulation in | 2 and 
phase. 

Now we turn to the numerical evidence for the persis- 
tence of DBs in the thermalized /3-FPU chain. As ob- 
served before [8], there are DBs present in the transient 
state. As nonlinearity increases, the duration of transient 
becomes shorter. After the energy redistributes among 
all the modes to achieve thermal equilibration, our simu- 
lations show that the spatially localized, high frequency 
excitations still exist. These DBs can interact with each 
other and may be destroyed by collision processes with 
other DBs or with the renormalizcd waves. The spatial 
structure of these excitations very much resembles the 
idealized breather oscillations in the absence of spatially 
extended waves: they "live" above the high frequency 
edge of the dispersion band and their lifetime is suffi- 
ciently long (on the order of 10-100 DB oscillations) to 
behave like a quasiparticlc. Note that, under certain con- 
ditions, supersonic solitons may arise from the /3-FPU 
system as another kind of localized excitations [15,16]. 
However, they were not observed in our thermalized sys- 
tem. Fig. 4 is the energy density plot which shows the 
time evolution of energy of each particle for the tran- 
sient (recording starting time T — 5x 10 2 ) and thermal- 
ized (T = 5 x 10 5 ) states, respectively. Fig. 4(c) and 
(d) display the energy as a function of site at To corre- 
sponding to Fig. 4(a) and (b) respectively. In the tran- 
sient case (Fig. 4(a)) the spatially localized objects (dark 
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stripes) that carry sufficiently large amount of energy are 
clearly observed. Fig. 4(c) is a snapshot of the energy 
density plot (4(a)) at To. Here the DBs are seen as lo- 
calized peaks [8]. After thermalization the spatial struc- 
ture looks different (Fig. 4(b)). The system now consists 
of the renormalized waves (straight cross-hatch traces in 
Fig. 4(b)). On the top of these waves, the localized struc- 
tures similar to DBs manifest themselves as the wavy 
dark trajectories (in Fig. 4(b)). Although the snapshot 
(Fig. 4(d)) of the energy density plot (Fig. 4(b)) indicates 
that in thermal equilibrium the energy is more evenly dis- 
tributed among particles, spatially localized structures 
are clearly observed. 

Since there are renormalized waves in the system, 
which also carry energy, we need to find a way to dis- 
tinguish between these waves and DBs. We use a fre- 
quency filter that cuts out the lower side of the Fourier 
spectrum and leaves the high frequency part unmodified, 
i.e., f{g(n,t)) = ^RF^ 1 (H UJ (F(g))), where 3? denotes the 
real part, F is a time Fourier transform, H u eliminates 
all frequencies below Lu cu t and g(n, t) is a dynamical vari- 
able that is being filtered. By applying this filter to the 
displacement q n to obtain = f(q n ), we can show the 
existence of DBs even for strong nonlinearities, for ex- 
ample, P = 25. Fig. 5(a) shows a clear example of a 
DB excitation reconstructed using the filtered q^ with 
Wcut — 7. Fig. 5(b) shows a typical spatial profile of the 
DB taken from Fig. 5(a), which strongly resembles the 
idealized DB [17]. Finally in Fig. 6, we present the evi- 
dence that there is a turbulence of DBs, which chaotically 
ride on renormalized waves. The corresponding energy 
density distribution along with the distribution of fil- 
tered displacement in a zoomed region qf(t) is displayed 
in Fig. 6(a) and (b), respectively. After the lower modes 
from the displacement q n are filtered, one can clearly ob- 
serve that the remaining high frequency oscillations are 
spatially highly localized, with the same characteristics 
as an idealized breather. The detailed time dynamics of 
the DB shows the main characteristics of breathers: the 
values of q^ change signs periodically (as indicated by the 
alternating white and black spots along the trajectory) 
as the DB moves in space, with a spatial span of 2 or 3 
sites only, as seen in Fig. 6(b). 

In conclusions, we have presented an interesting dy- 
namical scenario of the /3-FPU chains in thermal equi- 
librium: (i) For strong nonlinearity the linear disper- 
sion relation is effectively renormalized, which allows one 
to treat even strongly nonlinear systems as if they were 
weakly nonlinear; (ii) On top of renormalized waves, the 
strongly nonlinear system is also characterized by the 
turbulence of discrete breathers. 

Acknowledgments we thank Sergei Nazarenko for 
discussions. 




FIG. 1. Power spectra for the thermalized state for relatively 
strong nonlinearity (3 = l(dots), /3 = 8(circles) and f3 = 32(pluses). 
The thick solid line is for comparison. Note that time window 

T = 10 was used for averaging the power spectrum. 
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FIG. 2. (a) Density plot of spectrum for j3 = 1 and H = 200 
(ln(|afe(uj)| 2 ) is plotted), (b) renormalized linear frequency as a 
function of j3: analytical predictions (pluses) and numerical mea- 
surements (circles). /3 0,2 is shown for comparison (dashed), (c) 
effective nonlinearity as a function of /3 before (diamonds) and af- 
ter (crosses) renormalization (N=128). 
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FIG. 3. (a) and (b) |aj.| and phase respectively, using the bare 
linear dispersion ojj. = 2sin(7i-fc/7V). (c) and (d) |a&| and phase 
respectively using the renormalized dispersion ojk (Eq.3). Dashed: 
k = 1, solid: k = 20. The dispersion time scale is described by 
Tfc = 27r/ojfe for (a) and (b) and = for (c) and (d). (Note 

that according to our definition = |afc|e I ^+" fct ', the phase does 
not include the linear rotation.) 
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FIG. 4. Energy density evolution of the transient state (a) and 
(c), of the thermalized state (b) and (d), respectively ({3 = 1 and 
H = 200). In (a) and (b) the darker strips correspond to high en- 
ergy localizations, (c) and (d) are the snapshots of energy density. 
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FIG. 5. (a) Evolution of a discrete breather in thermal equilib- 
rium, (b) typical snapshot of the breather. /3 = 25, H = 200. 
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FIG. 6. Turbulence of discrete breathers (N=1024): (a) evolu- 
tion of energy density, (b) zoomed in q„ (t ) of the area indicated by 
the rectangle in (a). 
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